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Abstract. It has been shown that a Hamiltonian with an unbroken VI symmetry 
qq ' also possesses a hidden symmetry that is represented by the linear operator C. This 

symmetry operator C guarantees that the Hamiltonian acts on a Hilbert space with an 
inner product that is both positive definite and conserved in time, thereby ensuring 
that the Hamiltonian can be used to define a unitary theory of quantum mechanics. In 
£C} , this paper it is shown how to construct the operator C for the "PT-symmetric square 

(N ■ well using perturbative techniques. 
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£ . 1. Introduction 

cr 

The discovery PJ EJ E! that there were huge classes of PT-symmetric non-Hermitian 
Hamiltonians of the form H = p 2 + x 2 (ix) € (e > 0) whose spectra were real and positive 
led to the investigation of many new kinds of PT-symmetric model Hamiltonians. One 
particularly elegant model is the PT-symmetric square well, whose Hamiltonian on the 
domain < x < n is given by 

H = p 2 + V(x), (1) 

where V(x) = oo for x < and x > n and 

v(x) = l ie for I <a;<7r ' (2 ) 
V[X) -\-te forO<x<|. {Z) 

This Hamiltonian reduces to the conventional Hermitian square well in the limit as 
e — > 0. For H in the parity operator P performs a reflection about x = |: 
P : x —>■ TV — x. The PT-symmetric square- well Hamiltonian was invented and first 
examined by Znojil jl] and it has been heavily studied by many other researchers 

0E10IB]. 

The principal challenge in understanding non-Hermitian PT-symmetric Hamilto- 
nians was to show that they describe unitary time evolution. This was accomplished by 
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the discovery of a hidden symmetry operator called C. This operator is used to define 
the Hilbert-space inner product with respect to which the Hamiltonian is self-adjoint 
In Ref. |5] the C operator in coordinate space was shown to have a representation 
as a sum over the eigenf unctions (j> n {x) of the Hamiltonian: 

oo 

c ( x , V) = Y1 ^{x)(j) n (y), (3) 

n=0 

where the eigenfunctions are normalized so that they are eigenstates of the VT operator 
with eigenvalue 1, 

VT(f) n (x) =4> n {x), (4) 
and the integral of the square of the nth eigenfunction oscillates in sign: 

dx^x)} 2 = (-!)". (5) 



The discovery of the C operator led immediately to attempts to calculate it 
for various model Hamiltonians. For the elementary PT-symmetric non-Hermitian 
Hamiltonian H = |p 2 + ^x 2 + ix, the exact C operator is given by JU] 

C = e- 2p V. (6) 

However, for more complicated Hamiltonians the C operator cannot be obtained in 
closed form. It was shown in Ref. |H] how to use perturbative methods to evaluate the 
sum in (J3J) for the Hamiltonian 

H = p 2 + x 2 + iex 3 . (7) 

In Ref. jT2j this perturbative procedure was extended to quantum-mechanical 
Hamiltonians having several degrees of freedom. 

The perturbative methods used in Refs. JT] and ^2] were not powerful enough to 
be used in quantum field theory, so a simple recipe for finding C was devised that can 
be used in systems having an infinite number of degrees of freedom [TJIj. The procedure 
was to solve the three simultaneous algebraic equations satisfied by C: 

C 2 = l, [C,VT]=0, [C,H]=Q. (8) 

This recipe gives the C operator as a product of the exponential of an antisymmetric 
Hermitian operator Q and the parity operator V: 

C = e Q V. (9) 

Note that Q = —2p for the C operator in (jHJ). Mostafazadeh has shown that the square 
root of the positive operator can be used to construct a similarity transformation that 
converts a non-Hermitian "PT-symmetric Hamiltonian H to an equivalent Hermitian 
Hamiltonian h PJ: h = e~ Q / 2 He Q / 2 . 

In all the examples studied so far the C operator is a combination of integer powers 
of x and integer numbers of derivatives multiplying the parity operator V. Hence, the 
Q operator is a polynomial in the operators x and p = —i4z- The novelty of the PT- 
symmetric square-well Hamiltonian (0121) is that C contains integrals of V and thus the 
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Q operator, while it is a simple function, is not a polynomial in x and p and therefore 
cannot be found easily by the algebraic perturbative methods that were introduced in 
Ref. PHI- Thus, in Sec. El we calculate C for this Hamiltonian by using the perturbative 
techniques that were devised in Ref. jTT] . In Sec. H3 we make some concluding remarks. 



2. Perturbative calculation of the C operator 



The procedure we use here is as follows: First, we solve the Schrodinger equation 

-(f>'^x)+V(x)(P n (x) = E n (f> n (x) (n = 0, 1, 2, 3, . . .) (10) 

subject to the boundary conditions 0„(O) = n (7r) = 0. We obtain the eigenfunction 
4> n (x) as a perturbation series to second order in powers of e. The eigenfunctions are 
then normalized according to (|IJ) and (J2J). Next, we substitute the eigenfunctions into 
the formula Q and evaluate the sum. The advantage of the domain of the square well 
being < x < n is that this sum reduces to a set of Fourier sine and cosine series that 
can be evaluated in closed form. After evaluating the sum, it is convenient to translate 
the domain of the square well to the more symmetric region — ~ < x < |. On this 
domain the parity operator in coordinate space is V(x, y) = 5(x + y). Finally, we show 
that the C operator to order e 2 has the form in @, and we evaluate the function Q to 
order e 2 . Our final result for Q(x, y) on the domain 

Q(x, y) = \ie[x -y + e(x-y)(\x + y 

where e(x) is the standard step function 



2 < *~~~~ X *\ 2 IS 



TV 



)] + 0{e* 



'.ID 



six 



1 


(x > 0), 





O = o), 


-1 


(x < 0). 



(12) 



2.1. Solution of the Schrodinger Equation 

We begin our analysis by solving the Schrodinger equation (jlOj) in the right (x > |) and 
left (x < f ) regions of the square well: 

:|(l-(-l) n ) s [ n ( n _|_ l) x 



+ 



VRW = a n I « 2 



i,, , , ,„•, /7T x\ [ — l) n cos(n + l)x 1. . . sinm + l)x 



£ 2 (i+(-in i_z 

2 2 



+ G(e 3 )} 



(n + 1 

l(i + (-in f---^ cos ^+i)% 



4 



(n+ If 



2(n + l) 2 
x 2 7rx 7r 2 \ sin(n + l)x 
T ~ T + 16 J (n + 1) 2 



2 /' 



(13) 



+ 



= a n <! i2 (1 ( 1)n) sin(n + l)x 



,i (1+( _ 1) n ) x (-l) n cos(n + l)x j sin(riH-l)x 
2 (^TT) + 2( " { } } 2(n + l) 2 
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f(M-i) n ) 



1 xcos(n + l)x fx 2 7T 2 \ sin(n + l)x 

T + ( ~ ' 4 (n + 1) 3 + \ 8 ~ 16/ (n + 1) 2 



+ 0(e 3 



(14) 



These eigenfunctions and their first derivatives are continuous at x = ~. 

Having found the eigenfunction (fi n (x) to second order in e, we can give the formula 
for the corresponding eigenvalues: 



(n + 1 



(_l)n [2 _ ( _ 1 ) r j 2 ^ ( 

A(n+1) 2 



l -e 2 + 0(e i ). 



However, these eigenvalues are not needed for calculating the C operator. 



2.2. Normalization of the Eigenfunctions 

The normalization requirements in (|3J|HJ) give the value of the coefficient a n in (|13H14j) : 



l-(-l)' 



(2 



■l) n 7V 2 



e 2 + 0{e 4 



(15) 



(6 - 2(-l) n )(n + l) 4 16(71 + 1)2, 
With this normalization, the VT inner product between 4> m (x) and 4> n (x) is (— l) n 5 mn + 

2.3. Calculation ofC(x,y) to Leading Order (Zeroth Order) in e 

The next step is to construct the operator C(x,y), which is given in (jSJ) as a sum, by 
directly substituting the eigenfunctions (j> n {x) from (jlHH14j) . In general, there are four 
different regions of x and y to consider: 

(i) x > f, y > f, 



2 ' 

— 

2 ' 



x > 5, y < f : 



(iii) i<f,y<i : 

(iv) x < |, y > f . 

However, to zeroth-order in e, ra is common to all four regions and the calculation is 
easy. We find that 



*i 

C(°\x,y) = -y^(-l) n sin(n + l)x sin(n + l)y. 

71 / — ' 



(16) 



n=0 



This is just the Fourier sine series for the parity operator in the range < x < n: 

C^(x,y) = 8(x + y-n). (17) 

On the symmetric domain — | < x < | this formula becomes 

C(°)(x,j/) = 5(x + y), (18) 

which is equivalent to the coordinate-space condition of completeness. 
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2. 4- Calculation ofC(x,y) to First Order in e 



The calculation of C(x, y) to first order in e requires the evaluation of Fourier sine 
and cosine series. These are expressed in terms of single and double integrals of delta 
functions. Here, we describe the calculation for the region x > |, y > |. The calculation 
for the other three regions is similar. 

From (j!3|) and (fT3j) the first-order contribution to C(x, y) is 



C«(x,y) = i£ 

TT z — ' 



n=0 



(tt — x)- 



^ TT 



TT 

n=0 
l-(-l) 



7T 



(tt - y) 

oo 

E 



n + 1 
i(-l) 



cos[(n + l)x] sin[(n + l)y] 



ra=0 



— sin[(n + l)aj] cos[(n + l)y] 
sin[(n + l)x] sin[(n + 1 



{n + iy 



(19) 



The first two terms of can be expressed as single integrals of the parity operator, 
with the first having the upper limit x, 



dtS(t + y -tt) 



uo 

--E 

n=0 

- oo 



n=0 



71+1 
(-1)" 

n+ 1 



cos[(t2, + l)x] sin[(n + l)y] 
sin[(2n + 2)y] 



(20) 



and the second having the upper limit y, 



2 . \ 

dts(x + 1 — 7r) = — 



n=0 

oo 



n=0 



(-l) n 

— -j~Y sin[(n + l)x] cos[(n + l)y] 
n + 



y sin[(2n + 2)x) 



(21) 



The third term of involves sines of even x and y. To obtain a series of sines of 
even x and y, we subtract the series represented by the parity operator 5(x + y — tt) from 
the series represented by 5(x — y). The factor (n + l) 2 in the denominator implies that 
the third term of may be expressed as a double integral. To maintain the symmetry, 
we note that a double integral with respect to x contributes half of the third term of 

while a double integral with respect to y contributes the other half. 

This analysis allows us to express C^ 1 ' in the region x > |, y > | in terms of single 
integrals of the parity operator and double integrals in x and y: 



{ (f - i) />« + *-" 







pX 


+ 






4 


J-k/2 




X 


TT 


+ 






4 


' 4 



dt 



ds 



S{s-y) -8(s + y-Tr) 



w/2 

+ (x y) 



(22) 
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We simplify this result by evaluating integrals over delta functions and obtain 

C«(x, y) = \i{\x - y\ + x + y - 2ix) (x > f , y > f ) . (23) 

The calculation of C^(x, y) for the remaining three regions follows a similar 
procedure and we get 

C (1 \x,y) = 

\i[(x — 7r) 6{x + y — 7r) + y 9(tt — x — y)] 



\i[-\x-y\ + (x + y)] 
\i[{y — 7r) 6{x + y — 7r) + x 9(ir — x — y)\ 

where 9(x) is the Heaviside step function, 

1 (x > 0), 



9{x) 







(x > |, y < f) , 
(a; < f , y < |) , (24) 

l x < I > V > I ) > 



(25) 



(a; < 0). 

Finally, we condense the four expressions for C^\x, y) in the four different regions into 
a single expression: 

C {1) (x,y) = \i[x + y - n - 9(ir - x - y) (\x - y\ - tt) 
+ 9(x + y-n) (\x-y\ - tt)]. 

On the symmetric region — | < (x,y) < |, this expression becomes 

C (1 \x, y) = \i[x + y + e(x + y) (\x - y\ - tt)]. 

We plot the imaginary part of in Fig. ^ as a function of x and y. 



(26) 
(27) 



£.5. Calculation ofC(x,y) to Second Order in e 

The procedure for calculating C^(x,y) is similar to that used for calculating C^'(x,y), 
albeit more tedious. We must calculate sums of products of sines and cosines, but 
this time the presence of factors of (n + l) 4 , (n + l) 3 and (n + l) 2 in the denominator 
requires the use of quadruple, triple, and double integrals of delta functions to simplify 
the expression for C^(x,y). We discuss the calculation of C^ 2 \x,y) explicitly for the 
region x < ~, y < ~. The calculation of C^ 2 '(x, y) for the other three regions is similar. 
From (jEt and (fTKj) we see that the second-order calculation of C(x,y) gives 

x cos[(2n + l)x] sin[(2n + l)y] 
4 



C^(x,y) = -Y: 

7T z — ' 



n=0 

2 ^ 
+ - 



TT 



(2n + l) 3 
y sin[(2ra + l)x] cos[(2n + l)y] 



n=0 L 

oo 



+ -E 

7T Z — ' 

n=0 

r, OO 

--E 



y 2 



(2n + l) 3 
(— l) n sin[(n + l)x] sin[(n + l)y] 



n=0 



(n + 1) 2 

x cos[(2n + 2)x] sin[(2n + 2)y] 
4 (2n + 2) 3 
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Figure 1. Three-dimensional plot of the imaginary part of C^\x, y), the first-order 
perturbative contribution in (|27|l to the C operator in coordinate space. The plot is on 
the symmetric square domain — ^ < (x,y) < -|. Note that C^\x,y) vanishes on the 
boundary of this square domain because the eigenfunctions 4> n {x) in ||2Jl are required 
to vanish at x — and x = u. 



7T ^ 



n=0 
oo r 



--E 



n=0 
oo r 



--E 

7T Z -' 



n=0 
oo r 



--E 

n=0 



y sin[(2n + 2)x] cos[(2n + 2)y] 
4 (2n + 2) 3 

1 sin[(2n + l)x] sin[(2n + l)y\ 

2 (2n + l) 4 _ 

1 sin[(2n + 2)x] sin[(2n + 2)y\ 

2 (2ra + 2) 4 _ 

(— l) n cos[(n + l)x] cos[(n + 1) 

1 



n + 1 



(28) 



We have been able to evaluate each of these Fourier series exactly and to express the 
result as multiple integrals over delta functions: 

<• tt/2 



7.1' 



ptt/2 ft/ 

C {2) (x,y) = - \x 2 I dt I ds5(s + y-n) 

J X Jt 

tt/2 



dt 



tt/2 



ds 



tt/2 



dr 5(r + y — n) 



tt/2 



dt 



tt/2 



ds 



tt/2 



tt/2 



dr / dp8(p — y) 
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pn/2 ptt/2 

— \xy I dt I ds 5(s + t — n) + \x 2 y 

J y J x 

- ±xy7i + ±x 3 + {x++y). (29) 
Evaluating these integrals gives C^^x, y) for the region x < |, y < |: 

C {2 \x, y) = -±\x - y\ 3 + \x 2 y + ±x 3 - ±xiry + (x <-+ y). (30) 

The calculation of C^ 2 \x,y) for the remaining three regions follows a similar 
procedure. Combining the contributions from the four regions and transforming to 
the symmetric domain — | < (x,y) < |, we obtain the single expression: 

C {2 \x,y) = ^7r 3 + \xyii - j^n 2 (x + y)e(x + y) + |yr(a;|x| + y\y\) e(x + y) 

- + y 3 ) £ ( x + y)~ ^(y 3 - x3 ) £ (y - x ) 

- \xy{\x\[6{x - y) 9(-x - y) + 9(y - x) 9(x + y)) 

+ \y\[9(y - x) 9(-x - y) + 9(x - y) 9{x + y)]}. (31) 
We have plotted the function C^ 2 \x,y) on the symmetric domain — | < (x,y) < | in 

Fig.m 




Figure 2. Three-dimensional plot of C^ 2 '(x, y) in l|31|l on the symmetric square domain 
— 5 < (x,y) < ^. The function C^(x,y) vanishes on the boundary of this square 
domain because the eigenfunctions from which it was constructed vanish at the 

boundaries of the square well. 

In summary, our final result for the C operator to order e 2 on the symmetric domain 
— | < (x, y) < | is given by 

C(x, y) = 5(x + y) + eC« (x, y) + e 2 C^ (x, y) + 0(e 3 ) , (32) 
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where C^(x, y) is given in (|27jl and C^{x,y) is given in (j31j). We have verified by 
explicit calculation that to order e 2 this C operator obeys the algebraic equations (0)- 
For example, in coordinate space the third of these equations, C 2 — 1, reads 

/ ' dyC(x, y) C(y, z) = 5(x - z) + 0(e 3 ). (33) 

J-w/2 

2. 6. Calculation of the Q operator 

The last step in the calculation is to determine the operator Q from the result in (j3"2jl by 
using (JHJ). This is a long and difficult calculation: We first multiply C(x, z) on the right 
by 5(z + y), the parity operator V in coordinate space, and then integrate with respect 
to z. This gives the coordinate space representation of CV = e Q . Next, we take the 
logarithm of the resulting expression and expand it as a series in powers of e to obtain 
Q. We find that the coefficient of e 2 in this expansion is zero, and thus we obtain the 
simple result in (jlljl . which is the principal result in this paper. 



3. Conclusions 



In this paper we have used perturbative methods to calculate the C operator to second 
order in powers of e for the complex PT-symmetric square- well potential (j2j). Expressing 
C in the form e®V, we have found that the operator Q for this model has an expansion 
in odd powers of e, just as in the case of the cubic "PT-symmetric oscillator whose 
Hamiltonian is given in (|7|) . Our result (jlljl for Q is an elementary function. We 
have verified our calculation of the C operator by showing that it satisfies the algebraic 
conditions (jHJ). 

The most noteworthy property of the C operator is that the associated operator Q 
is a nonpolynomial function, and this kind of structure had not been seen in previous 
studies of C. At the beginning of this calculation we expected that for such a simple 
PT-symmetric Hamiltonian it would be possible to calculate the C operator exactly and 
in closed form. We find it surprising that even for this elementary model the C operator 
is so nontrivial. 
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